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Abstract 

We discuss optical potentials for the nuclear elastic scattering from a variational viewpoint. 
Density-dependence in the effective N-N interaction leads to density rearrangement terms, in 
addition to the conventional folding term. Effects of the rearrangement on the N-A optical potential 
are illustrated in the nuclear-matter limit. Closely relating to consistency with the saturation, 
the rearrangement appreciably improves the isoscalar optical potential depth over the previous 
folding model calculations. The rearrangement gives stronger effects as the density grows. We 
also present rearrangement terms in the A-A double-folding potential. Since the rearrangement 
terms are relevant to the nuclear structure but should be handled within the reaction model, N-N 
effective interactions applicable both to structure models and to the folding model will be desired 
for unambiguous description of the nuclear elastic scattering. 

PACS numbers: 24.10. Cn,24.10.Ht,25.40.Cm,25.70.Bc 
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The optical potential for elastic scattering is basic to describing various nuclear reac- 
tions 1]; the direct reactions including inelastic scattering, charge exchange and nucleon 
transfer reactions, which are treated in the distorted- wave Born approximation (DWBA), 
as well as the compound nuclear reactions. The folding model has successfully been ap- 
plied to obtain the optical potentials, their real parts in particular, based on effective N-N 
interactions 0. 

In regard to the effective N-N interaction, density- dependence is required to reproduce 
the saturation of nuclear densities and energies, as qualitatively disclosed by the Brueck- 
ner theory [31]. However, Brueckner's g-matrix does not reproduce the nuclear saturation 
quantitatively (as long as we apply it with two-body interaction within the non-relativistic 
framework). In most nuclear structure problems we use phenomenological effective inter- 
actions in which density- dependence is contained so as to reproduce the saturation j^J. 
Density-dependence has also been introduced in the effective N-N interaction for the fold- 
ing model 5]. Efforts have been made to connect the folding- model interaction to the 
saturation 6, Q]. It seems hopeful to unify the folding model at relatively low incident 
energies and nuclear structure models such as the mean- field approximation ^, 0|. How- 
ever, there exists inconsistency between them in treating the density- dependence. In this 
article we reconsider the folding model from a variational viewpoint, and show where the 
inconsistency comes from and how it should be resolved. 

We here express an effective N-N interaction in the following form: 

V = J2 C n[p]w n , (1) 

n 

where w n is a two-body interaction operator, C n [p] its density-dependent coupling constant. 
As in Refs. , the coupling constants are assumed to depend on the total density p. 

Although there could also be energy- dependence in the coupling constants, it does not play 
an essential role in the arguments below. 

We first consider the N-A scattering. The system has (A + 1) nucleons, with the total 
wave-function |^a+i) under the total Hamiltonian H A +i- We decompose I^a+i) as 

I^A+l) ~ |*A> ® VI>n) , (2) 

where \^f A ) is the wave-function of the target nucleus A, and \4>n) denotes the wave-function 
of the incident and scattered nucleon. In the folding model, we customarily assume \^a) to 
be the ground-state function l^^), which satisfies HaI^a) = /• This assumption 

is practically equivalent to the frozen density approximation. With the decomposition (j2J), 
the Schrodinger equation for ip^ (or the Dirac equation in relativistic approaches, with 
appropriate subtraction of the mass) is derived by the variational equation, 



8ip 



(* A+1 \(H A+1 -E)\* A+1 ) =0, (3) 



N 
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which can be written as 



h N \ip N ) = En I ipN ) ; h N = K N + U, E N = E-E A , (4) 

where Ea = (^a\Ha\^a) (~ E A ^), stands for the kinetic energy operator and U is 
regarded as the optical potential. In applying Eq. (J2J) to (jSJ), one should notice that pa+i 
obtained from |^a+i) should be used in C n [p\ of Ha+i, while Ha contains pa which is 
calculated from ). Having certain contribution to U, 5p = Pa+i — Pa = ^n^n is n °t 
negligible though infinitesimal at each position. We then obtain 

U\M = £ \c n \p A ] • (* A |to n (|* A > ® \M) + S( r -^ - r N ) • |^)| . (5) 

Here (O)a = (^a\0\^a), and 5{^p- — r N ) equates the center-of-mass position of the 
interacting nucleons in A to the position of the scattered nucleon. The first term on the 
rhs in Eq. (|SJ) is the usual folding potential. It is remarked that the density rearrangement 
term (DRT) comes out as the second term. While such rearrangement terms are known in 
structure models and in some reaction approaches (e.g. the one within the time- dependent 
Hartree-Fock approximation), the DRT has not been taken into account in the folding model. 
We call U in Eq. © renormalized folding potential (RFP), because of the correspondence to 
the renormalized Brueckner-Hartree-Fock approach 0| as discussed later. The RFP differs 
from the conventional folding potential because of the presence of the DRT. 

To illustrate significance of the DRT, we consider the nucleon scattering on the uniform 
nuclear matter with equal numbers of protons and neutrons. This nuclear-matter folding 
potential corresponds to depth of the isoscalar optical potential in the N-A scattering. As 
far as contribution of the dynamical polarization is negligibly small, the RFP is 

U(k; k F ) = V \c n [p\ -^-r V / d 3 k' (karkVVK|karkW) + ^ ^M^)] , 



dp 



(6) 

where k-p denotes the Fermi momentum, and Q is the volume of the nuclear matter. We 
have p = (2/37r 2 )fcp, as usual. Note that U becomes independent of the spin and isospin 
component a and r in the symmetric nuclear matter, as well as of direction of the momentum 
k. Obviously U(k;k F ) is connected to the single-particle (s.p.) energy. The s.p. energy is 
defined by variation of the nuclear matter energy E with respect to the occupation number 
n kar, which is unity for k < kp and vanishes for k > kp. The nuclear matter energy is 
expressed by 

E = -, — tt7 / o9k n\r aT 

1 Q 2 r 

+ 2^^(2T ^ J d 3 kd 3 k' (kaT'k'a'T'\w n \kaT'k'a'T , )n k(TT n k , a/T , , (7) 
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and therefore 

5E k 2 

This e(k; hp) is equivalent to in Eq. (0J) in the scattering problems, manifesting similarity 
of the folding model to the mean-field picture. Let us now recall the Hugenholtz-van Hove 
(HvH) theorem The saturation point is obtained as the minimum of £ = E/A = E/pQ. 
Since the Fermi energy is given by e F = e(k F ; k F ) = (l/fl) (dE /dp), we obtain 

SF = £ + P dp~- (9) 

This yields e F o = £(&fo; &fo) = £mm at the saturation point, where k F0 represents the Fermi 
momentum at the saturation and £ m j n the saturation energy. We then have 0| 

U(k Fo] k FO )=£ min -^. (10) 

Since the saturation is directly linked to the variation of £, this relation is essential to 
consistency with the saturation. Furthermore, the effective mass at the saturation point, 
Mq, is given by 



kpo = de(k; k F0 ) 
M * ~ dk 



k F0 dll(k; k F0 ) 



M dk 



(11) 

fc=fcpo 



which derives 



dU(k; k F0 ) 



dE 



N 



M* 

= 1-^2.. (12) 

t * M 

Here E N = e(k;k F ) is related to k via Eq. (JSJ). Note that, even without energy-dependent 
coupling constants, non-locality in the interaction gives rise to /c-dependence of U, leading 
to -Eat- dependence. For the scattering problems, E N must be positive. However, with 
the extrapolation k — > k F0 (i.e. E N — > £ min m — 16MeV), U(k;k F0 ) is constrained by 
Eqs. ()10pi2jh as long as it is continuous. 

Figure H] shows U(k; k F ) of Eq. (jHJ) as a function of En at p = 0.16 fm~ 3 . The calculations 
are implemented with several interactions used in the mean-field approaches as well as 
with a density-dependent interaction developed for the folding model, BDM3Y1 jo]. In the 
mean-field interactions, a p-dependent coupling constant is introduced for a contact term 
of the interaction. All of their results meet almost at the same point by extrapolation to 
E^ ~ — 16MeV, fulfilling the HvH theorem. They also give slopes close to one another at 
E N ps — 16MeV, corresponding to Mq ps 0.7M jslE3|- However, as E N increases, U(k; k F ) 
behaves differently among various effective interactions. This is not surprising because these 
mean-field interactions are determined from the structural information, which is not sensitive 
to U(k; k F ) far from the saturation point. The Skyrme interaction 0,^^ necessarily gives 
linear dependence on E N . While D1S [15] does not give reasonable E N - dependence, Dl [l^ 
and M3Y-P2 |9| have E N - dependence similar to that of the empirical depth of the optical 
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FIG. 1: Energy-dependence of the folding potential for the nucleon scattering on the nuclear 
matter, at p = 0.16 fm -3 . Results from several N-N effective interactions are compared; Skyrme 
SLy4 Q (green line) Gogny Dl 16] (blue dashed line) and D1S 15] (blue solid line), M3Y-P2 
(red line), BDM3Y1 |6| without the DRT (black solid line) and that with the DRT (black dotted 
line). Experimental values are represented by squares, which are obtained from depth of the real 
part of the isoscalar optical potential 
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potential. There is a class of effective interactions in which the M3Y interaction 
multiplied by a ^dependent factor 00,0. We here treat BDM3Y1 a 8 a representee of 
them. For BDM3Y1, both the results with and without the DRT are presented. We do not 
use the factor giving additional energy- dependence (denoted by g(E N ) in Ref. [(J). Although 
BDM3Y1 is adjusted so as to reproduce fcpo and £ m i n , the folding potential calculated with 
BDM3Y1 in Ref. is contradictory to Eq. dTUJ) , because the DRT is discarded. If the DRT 
is added, depth of the folding potential with BDM3Y1 is much closer to the empirical values, 
and is almost indistinguishable from that with M3Y-P2. Thus the RFP improves the optical 
potential depth over the previous folding model approach to a substantial extent, due to 
the presence of the DRT. Even the additional -E^-dependent factor seems unnecessary for 
BDM3Y1 at En < 200 MeV. This result demonstrates significance of the constraint fllOl) . 

In Fig. |2l p-dependence of U{k;k^) is depicted for several En. Because of the in- 
dependence presented in Fig^ the folding potential becomes less attractive, or even repul- 
sive, as En increases. As viewed in the BDM3Y1 results, effects of the rearrangement are 
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FIG. 2: Density-dependence of the folding potential for the nucleon scattering on the nuclear 
matter at En = 0, 60, 120 and 180 MeV (from the upper line to the lower), a) Results from 
BDM3Y1, with (dotted line) and without (solid line) the DRT. b) Results from M3Y-P2, for which 
the DRT is taken into account. 



the stronger for the higher p. Therefore, although the rearrangement will not seriously in- 
fluence highly peripheral reactions, it affects the volume integral of the potential. It will be 
important to reassess reaction calculations by applying the RFP. Whereas BDM3Y1 with- 
out the DRT gives attractive U(k; kp) even at twice the saturation density, the RFP comes 
repulsive at high p for relatively high E^. This p-dependence is qualitatively similar to that 
with M3Y-P2. Another interesting point suggested by Fig. El will be that the DRT tends 
to produce or enhance the so-called wine-bottle-bottom shape in the optical potential at a 
certain energy region, which significantly influences the analyzing power 19 1. 

Whereas we have assumed the interaction form of Eq. ((TJ, similar arguments hold for the 
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FIG. 3: Goldstone diagrams regarding the optical potential. Left: lst-order diagram in the hole- 
line expansion. Middle: 2nd-order diagram due to rearrangement for the Pauli exclusion operator. 
Right: 2nd-order diagram due to rearrangement for the energy denominator. Exchange diagrams 
are not shown. Solid and dashed lines represent the nucleon and the g-matrix, respectively. 



microscopic N-N interaction obtained in terms of the g-matrix lfj, [12J, |2fJ, |2l|, |22j . From 
the g-matrix viewpoints, the density- dependence in the effective interaction arises from the 
Pauli principle, which affects the Pauli exclusion operator and the self-energy in the energy 
denominator. In the N-A scattering, the projectile shifts the energy of the target nucleus 
because of the Pauli blocking; the s.p. state occupied by the projectile is blocked in addition 
to those occupied by the nucleons in the target nucleus. This additional blocking leads 
to a DRT, which is equivalent to the DRT in Eq. The density rearrangement in the 
(/-matrix was discussed in Refs. jl2l l2o| . In terms of the hole-line expansion, the density 
rearrangement is primarily represented by the 2nd-order diagrams (rearran gem ent for the 
Pauli exclusion operator and for the energy denominator) shown in Fig. El |2fJ, while the 
conventional folding takes account only of the lst-order diagram. In this regard the RFP 
is compared to the renormalized Brueckner-Hartree-Fock approach in R ef. llCw] . Relevance 
of the rearrangement to the HvH theorem was also pointed out in Ref. [20|. However, its 
effects and importance in the optical potential have not been recognized sufficiently so far. 
The present results indicate that the DRT is not negligible. The DRT does not give all 
the 2nd order terms in the hole-line expansion. It seems that, rather than the order in the 
hole-line expansion, the constraint of Eq. plays an important role. 

We next turn to the A- A scattering. Let us denote the target nucleus and the projectile 
by A\ and A 2 , respectively. In this case we have the following DRTs in the RFP, 



= £ {({Cn[p Al + PA 2 ] ~ C n [ PAl _ 

n v 



w, 



C n [pA 1 + PA 2 ] ~ C n [pA 2 



w. 



(13) 
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The first term on the rhs in Eq. (JT3j) is evaluated only from | ) , with displacement of the 
p-dependent coupling constant in the interaction of Eq. ([TJ; vice versa for the second term. 
As the A 2 nucleus approaches Ai, the density changes, causing the energy shift through the 
p-dependence in the interaction. Note that, in the pa 2 ~^ limit, the first term recovers the 
DRT in Eq. ©. 

To discuss the density rearrangement in the A- A scattering from the ^-matrix viewpoints, 
we shall assume that the momentum is approximately a good quantum number in the s.p. 
states, for simplicity. If the A± nucleus is isolated, nucleons in A\ occupy the s.p. levels of 
k < ^f,Ai, and this distribution gives rise to the Pauli blocking effects among the nucleons 
in Ai. When the A 2 nucleus is present nearby, the s.p. levels of |k — K| < A;f,a 2 are a ls° 
blocked (K is the relative momentum of A 2 to A\). This additional blocking affects the 
A\ energy, producing a DRT which corresponds to the first term in Eq. (jl3j) . We note 
that the Pauli blocking due to the two Fermi spheres in the A- A collision was considered 
in Ref. 0, although relation to the folding model was not elanded. It is often a good 
approximation to replace the Pauli blocking effects by a function of the total density 
Under this approximation, we go back to the same arguments as those based on Eq. (JIJ), 
which have yielded the DRTs in Eqs. ()5I13|) . 

The DRTs represent energy shift of the nuclei during the scattering process. They should 
be handled within the reaction model, because they affect the nuclear properties dynam- 
ically. Now the following question is raised: should the interaction v in Eq. (JT3~|l be the 
interaction in the structure model or that in the reaction model? Conventionally, separabil- 
ity of nuclear reaction problems from nuclear structure has been postulated, and different 
effective interactions have been used between the nuclear reaction theory and the nuclear 
structure theory. However, there are many cases in which reaction problems strongly couple 
to nuclear structure. Unified description using consistent interactions is desirable for total 
understanding of the phenomena. The above question suggests that the DRTs are located at 
an intersection of the structure theory and the reaction theory, exemplifying importance of 
unified description. To avoid ambiguity, consistent effective interactions should be applied 
to structure and reaction calculations. If we restrict ourselves to reactions at relatively low 
energies, it might be possible to develop an interaction applicable both to structure and 
reaction calculations. Analysis like that in Fig. ^may be useful for first selection of effective 
interactions for that purpose. 

In low-energy nuclear reactions, p-dependence of the effective interaction is important, 
and could somewhat be correlated with the saturation. In contrast, p-dependence is not so 
significant in high-energy reactions. If the incident energy is sufficiently high, the impulse 
approximation works well. In such energy regime, the phenomenological effective interac- 
tions employed in Figs. Q and El will lose their validity. Furthermore, effects of the density 
rearrangement are expected to be weak. It will be interesting to estimate quantitatively by 
the g-matrix approach how large the DRTs are at varying E^. 

Many-body forces have similar effects to the density- dependent two-body interactions 0| . 
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When a many-body force is introduced instead of the density-dependent interaction, we have 
DRTs again. For instance, in the N-A scattering with a three-body force, the interaction 
among the projectile and two nucleons in the target nucleus gives rise to energy shift of the 
target nucleus in an effective manner, which produces a DRT as in Eq. (JHJ). 

Effective interactions for inelastic scattering will also be subject to the density rear- 
rangement effects. Because the density-derivative of the optical potential is relevant to 
the inelastic scattering to collective states within the DWBA framework 24| . terms up to 
S 2 C n [p]/S 2 p come into the issue of the N-A process. This coincides the residual interaction 
in the random-phase approximation J^. Details will be discussed in a future publication. 

The folding model using density- dependent interactions has been successful in describing 
nuclear reactions to a certain degree, without the DRTs. There could be several reasons for it. 
First, some reactions may be insensitive to the rearrangement effects, because, e.g., they are 
highly peripheral. Second, some of the density-dependent interactions that have been used 
in the folding model were fitted to the reaction data. This may mask some rearrangement 
effects. Furthermore, the dynamical polarization effects should be treated properly, for full 
understanding of the optical potential. Still it will be important to reinvestigate individual 
reactions by taking the DRTs into account. As stated already, this may open a way to 
unified description of the nuclear elastic scattering and the nuclear structure. 

In summary, we have discussed the N-A optical potential based on a variational equa- 
tion. Density-dependence in the effective N-N interaction gives a density rearrangement 
term, which has been discarded in the conventional folding-model calculations. We call 
the optical potential with the rearrangement term renormalized folding potential. If the 
rearrangement term is taken into account, depth of the N-A folding potential is in better 
agreement with the experimental data than in the conventional folding model approach. 
We emphasize significance of a constraint due to the Hugenholtz-van Hove theorem, which 
assures consistency with the saturation. Effects of the rearrangement become stronger as 
the density increases. Rearrangement terms in the A- A scattering are also discussed. The 
rearrangement terms seem to be located at an intersection of the nuclear structure mod- 
els and the folding model. It is desired to explore effective N-N interactions applicable to 
both of these models in a consistent manner, for unambiguous description of the low-energy 
nuclear elastic scattering. 
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